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ABSTRACT 

A review is presented of the past 165 years of observation of the 4-day Cepheid Polaris, includ- 
ing the exciting results of the last 50 years, an interval that has produced three orbital solutions 
for the spectroscopic binary subsystem, recently resolved by HST, parameters for the optical 
companion, precise measurement of the star's trigonometric parallax and angular diameter, ev- 
idence for a rapid increase in its pulsation period, and observations of the dramatic decline and 
recent partial recovery of its light amplitude. There has been considerable discussion about the 
exact nature of the star, with potential resolutions summarized here. It is also noted that many 
of the star's characteristics are shared by a small number of other Cepheids that display rapid 
period increases identical to those predicted for stars in the first crossing of the instability strip, 
small light amplitudes, and intrinsic colors typical of variables lying near the center of the strip, 
where Cepheids of largest amplitude reside. While all members of the group appear to display 
the canonical traits of first crossers of the instability strip, Polaris has one unique peculiarity: a 
brief hiatus in its monotonic period increment between 1963 and 1966 during which the pulsation 
period underwent a dramatic decrease. Has the average brightness of the Cepheid also increased 
over the last millennium? 

Subject headings: stars: binaries: general — stars: variables: Cepheids — stars: individual: Polaris 



1. Introduction 

Polaris (a UMi, spectral type = F7 lb) is the 
nearest and brightest of ~ 4 classical Cepheid s 
detectable without optical aid ( Turner et al.l2 005). 
and possibly the most enigmatic. It has a pulsa- 
tion period of 3 d .97 and the se cond smallest light 
amplitude (ITurner et al. Il2009bh. Its light variabil- 
ity was suspected by Seidell ( 18521 ) from measures 
he made using a Steinheil eyepiece photometer 
that were corrected for atmospheric extinction and 
had a precision of about ±0 m .05. The accuracy of 
measures made with th e photometer is somewhat 
uncertain ( Zissell 1998), and in any case awaited 
the full development of the magnitude scale. Vari- 
ability was also suspected in visual comparisons 
of the star with [3 UMi made b y Schmidt bet ween 
1843 and 1856 (1144 estimates. ISchmidtlfl857h and 
Pannekoek betwe en 1890 and 1894 (510 estimates, 
Pannekoeklll913l) . but confirmation of the star's 
light variations would await photographic obser- 



vations by iHertzprung ( 1911 ). after the star's 4- 
day cycle of radial ve locity variability had been 
announced previously ( Campbel]||l899| ). The orig- 
inal light amplitude was very small, only ~ 10%, 
which is why it remained unnoticed for so long. 
In the General Catalogue of Variable Stars it was 
listed as a Type II C epheid as recently as 1970 



Kukarkin et al. 1970). although its metallicity 



jusenko et al.ll2005l ) and distance (see below) con- 
firm it to be a classical Cepheid. 

Polaris i s a t riple system (jKamper 1996; 
Evans et al.1 [2008) consisting of an unresolved 



F6 V dwarf, Polaris Ab, in a 30-year spectro- 
sco pic binary orbit and a F3 V d warf (Polaris 
B) (jTurnerlll977t lEvans et al.ll2008l ). 18" distant, 
as a visual physical companion. Polaris has a 
preci se Hipparcos parallax of 7r a ,h s = 7.72 ± 0.12 
mas ( van Leeuwen et al. 20071 ) that yields a dis- 
tance of 130 ± 2 pc. The corresponding lumi- 
nosity, (My) = —3.62 ± 0.05, implies overtone 
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(OT) pulsation, consistent with the Cepheid's si- 
nusoidal light curve. But there is no indication 
in Fourier analyses of photometric or radial ve- 
locity measures of a signal at the p utative funda- 
mental mode (FM) p eriod of 5 d .7 (iTurnei) l2006t 
ISprecklev fc Stevens! 120081: IBruntt et alj (20081 ). 
and a sinusoidal light cu rve may no t always 
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imply overtone pulsation (jSandage et al. 
Classification through Fourier parameterization 
is not necessarily reliable for a light curve of ex- 
tremely small amplitude consistent with a pure 
sine wave, and a preliminary investigation of 
nine other si nusoidal Cepheids a ssociated with 
open clusters ([Turner et al. 2006a) indicates that 
only four are likely OT pulsators. V1726 Cyg in 
the cluster Platais 1 is a good example of what 
might be the case for Polaris: a low-amplitude, 
sinusoidal Cepheid pulsating in the FM with 
(M v ) = -3.08 ±0.07 (|Turner et al.ll2006bl ). What 
is the situation for Polaris? 

2. Background Information 



Polaris 



( Witt et al 



illuminates 
2008h 



as 



diffuse 
well 



as a 



dust nebulosity 
reflection neb- 
ula ~ 1° south (jZagurv et al J Il999h . presum- 
ably the associated dust clouds accounting for 
the Cepheid's small reddening of Eb-v — 0.02 ± 



0.01 dGautheir fc Fernidll97ajUsenko et al.ll2005t 
Turned l2006F ~ The star has been the object of 
numero us ground-based inv estigations of its vari- 
ability (jTurner et alj 12005). ranging from visual, 
photographic, and photoelectric b rightness esti- 
mates to radial velocity measures fa oemer Il965t 
Kamper 19961 : Turner et al.l l2005h . A selection 
of some of the best photometric observations is 
shown in Fig. [TJ along with a variety of differ- 
ent radial velocity measures. Not illustrated are 
the many low-quality observ ations that have be en 
made over the past century (jTurner et al.ll2005[ L 

Recent years have witnessed great improve- 
ments in the quality of both the photometric 
and radial velocity monitori ng of Polaris, includ- 
ing observations from sp ace ( Sprecklev k. Stevens! 



2008: IBruntt e t al. 2008) an d high resolution spec- 



troscopy (lKampejll996l: lHatzes fe Cochranl feoOO: 
Turner et al.ll2006cl : iLee et al.ll2008h . High quality 
radial velocity measurements are probably the op- 
timum method for updating information on the 
period and amplitude changes of Polaris from 



ground-based observatories, although it is from ex- 
actly such data that evidence has been presented 
for an additional periodicity in the variability of 
Polaris, other than that arising from its pulsa- 
tional and orbital motions. Is the 119-day pe- 
riodicity found by ILee et al" ( 2008! ) evidence for 
a fourth star in the system, or is it attributable 
to seasonal effects on ground-bas ed equipment? 
A sim ilar f eature was detected bv | Dinshaw et al 
(|l989t ) and lHatzes k Cochran! (|2000|) in their ra- 



dial velocity studies of Polaris, but with a shorter 
period near 40 days. 
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Fig. 1 . — Photometric light curves f or Polaris (left) 
from: IStebbins fc Whitfordl Jl938h . de Vaucouleursl 
dl947l), iFeltz fc McNamaral dl980l ). iHendenl l|l980l ). 
and lBrown fc Bochonkd (1 19941): pha s ed ra di al veloc- 
ity va riation s (right) from Roemerl |l965l) . iKampei! 
(| 19961 ). and iTurner et all l|2006d ). with source tele- 
scopes indicated. 



3. Period and Amplitude Changes 

A recent study of th e period changes in Po- 
laris (|Turner et all 12005! ) updated a nd expanded 
upon more limited investigations by Ferro ( 19831) 
and others dFernie et al.lll993l : [Brown fcB o chonko 
ll994HEvans et alj|2002h . The burner et all (|2005l ) 
study included light travel time corrections for the 
orbi tal motion of P olaris about its F6 V compan- 
ion ( Kamper 19961 ). although it was noted in the 
analysis that such corrections did not appear to 
remove the effects entirely. The original analysis 
was redone for the present study with the inclusion 
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of recent observations, but excluding the original 
light travel time corrections, in order to determine 
if they could be established independently. The 
results are presented in Fig. 

The long-term trend in the pulsation period of 
Polaris from O-C analysis is that it is constantly 
increasing, except for an unusual "glitch" circa 
1965, during a gap in the observational coverage 
between 1963 and 1966. Prior to the glitch, the 
rate of period increase is calculated to have been 
4.46 ± 0.03 s yr _1 . Following the glitch it is cal- 
culated to be 4.19 ± 0.13 s yr _1 . Both values are 
consistent with results from stellar evolutionary 
model calculations for a star in the first cr ossing 
of the instability strip (jTurner et al.ll2006dl ). pro- 
vided that Polaris is a FM pulsator. If the Cepheid 
is in another crossing of the instability strip, say a 
third or fifth crossing, the observed rates of period 
increase are several ti mes larger than pre dicted 
from model predictions (jTurner et al ]|2006dl) . The 
shift in the O-C data during the glitch corre- 
sponds to an abrupt decrease in pulsation pe- 
riod, which, according to simple calculations tied 
to the period-density relation, could be accom- 
modated through the sudden acquisition by the 
Cep heid of about seve n Jovian masses of mat- 
ter (jTurner et al.ll2005l) . The true explanation re- 
mains a mystery, although it is interesting to note 
that the assimilation of any pre-existing planetary 
companions to Polaris is most likely to occur dur- 
ing the first crossing of the instability strip, when 
the star is growing to red supergiant dimensions 
for the first time. 
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Fig. 2. — O-C variations of Polaris from photometric 
and radial velocity data, without correction for orbital 
motion. 

The amplitude variations in V of Polaris are 
plotted in Fig. [31 which includes radial velocity 



amplitudes scaled by the usual factor of 50 km 
s _1 magnitude -1 . The light amplitude prior to the 
1965 glitch appears to have been undergoing a slow 
but steady decline, wit h the much more rapid de- 
cline noted by others (lFerrolll983l: iDinshaw et all 
ll989l:lBrown fc Bochonkolll994l:lEvans et alfeooj ) 
beginning roughly a decade later, reaching mini- 
mum at AV = m .025 circa 1988. It is presently 
increasing, but may decay completely by about 
2400 if the pre-glitch decline continues. 
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Fig. 3. — F-band amplitude variations of Polaris 
from photometric (open circles) and radial velocity 
(filled circles) data, with uncertainties indicated. The 
straight line is a weighted fit to the pre- 1966 data. The 
lower portion expands the data for the past 50 years. 

A longer-term overall brightening of ~ l m .7 
(from 3.6 to 1.9 in V from ~ 100 CE to the 
present) has been suggested (lEngle fc Guinan 



120041 : ICroswelll[2006l) . based upon an analysis of 
archival and recent brightness estimates. Such 
a gradual brightening over the last millennium 
would make Polaris unique among Cepheids, but 
appears to conflict wit h two sets of observations: 
those of lSchmidtl (|l857l ) that imply (V) ~ 2.0 pro- 
vided there is no color term in his c omparison of 
Polar is with the mu ch redder (3 UMi (jTurner et al 



20051 ). and those of iTurner et all (|2009aD over the 
past seven years that imply (V) — 1.991 ± 0.004 
s.d., a value which includes slight contamination 
from the light of Polaris B amounting to less than 
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m .006 (i.e., (V) ~ 1.997 for the Cepheid alone). 
In other words, the mean brightness of Polaris ap- 
pears to have been constant at V ~ 2.00 at least 
over the past century and a half. 

4. Basic Parameters 



A model a tmosp here analysis of Polaris by 
(|2005h finds an abundance pattern 



Usenko et al.l 
with [C/H] = 
-0.00, and 



= -0.17, [N/H] = +0.42, [O/H] = 
[Na/H] = +0.09, with the nitrogen 
and sodium enhancement and carbon (and oxy- 
gen) depletion being a strong signa t ure of CNO- 

argue 



processed elements. Usenko et al. ( 20051 ) 



that Polaris is a post red supergiant dredge-up 
star because of such a signature. Yet the pres- 
ence of CNO-processed elements in Cepheid atmo- 
spheres appears to be unrelated to red supergiant 
dredge-up, which does not take place in some evo- 
lutionary models, the contaminatio n occurring in 



stead in late main-sequenc e phases (Macder 2001 
Turner fc Berdnikovll2004[ ). The most likely mech 
anism is meridional mixing in rapidly-rotating B- 
type stars. As noted earlier, the rapid rate of red- 
ward evolution indicated by its rate of period in- 
crease suggests that Polari s is in the first cross - 
ing of the instability strip (jTurner et al.ll2006dr ). 
The rate is almost two orders of magnitude larger 
than predicted by stellar evolutionary models for 
a third crossing, although a fifth crossing might be 
an outside possibility. 

The angular diameter of Polaris has been mea- 
sured with the Naval Prototype Optical Inter- 



ferometer (jNordgren et al. 2000), yielding best 



estimates of #ud = 3.14 ± 0.02 mas and 6>ld = 
3.28 ± 0.02 mas. Cepheids are known to obey a 
tight period-radius (PR) relation established from 
Baad e-Wesselink and surfa c e brightness meth- 
ods (|Lanev fc Stobiel 1 19951: iGieren et all Il998t 
Turner fc Burkell2002l) . with (R) ex P§ . Cepheids 
also appear to obey a period-mass (PM) rela- 



tion established from cluster membership (jTurner 
19961). with M o c Pi For the PR relation of 
Turner fc Burke ( 20021 ). the implied distance to 
Polaris from its angular diameter is 93 + 2 pc 
(6>ld) or 97 ± 2 pc (6»ud) for FM pulsation, and 
122 ± 3 pc (0 L d) or 128 ± 3 pc (6> UD ) for OT pul- 
sation. For comparison, the Hipparcos parallax 
(jvan Leeuwen et al.1 120071 ) corresponds to a dis- 
tance of 130 ± 2 pc, while association with Polaris 



B implies distances of 101 ± 3 pc dTurner 20061 ) or 
109.5 pc (|Usenko fc Klochkovall2008l ). The limb- 
darkened angular diameter of Polaris does not 
appear to produce a distance estimate via this 
prescription that is entirely consistent with either 
of the values established by its cluster or trigono- 
metric parallaxes, which is puzzling. 

Polaris has a measured surface gravity of 
logg = 2.0 ± 0.3 (jUsenko et al.ll2005l) from LTE 
model atmospheres. The corresponding mass es- 
timates for Polaris using the relations described 
above are given in TablelU All of the estimates ap- 
pear reasonable, and imply a mass for the Cepheid 
lying somewhere between 4 and 7 -Mq. 

5. The Polaris Cluster 

The field surrounding Polaris appears to coin- 
cide with the remains of a poorly-populated star 
cluster. For example, the radial velocities of Po- 



tical: —16 km s 1 


Polaris, Roemerll 19651 Kamper 


1996: Turner et al. 


2006d). -15 km s" 1 (Polaris 


B. Kamperlll996HUsenko & Klochkova 


2008), and 


-11 km s" 1 (HD 5914, Abt & Biggs 


Il972h. A 



color-magnitude diagram for possible cluster stars 
from Hipparcos photometry is shown in Fig. [U 
Most of the stars, other than Polaris, have par- 
allaxes corresponding to distances of ~ 100 pc. 
The age isochrone for t = 8 x 10 7 years shown in 
Fig. |4] was chosen to correspond to Polaris as a 
first-crossing Cepheid. 

Simila r results are obta ined using 2MASS pho- 
tometry ( Cutri et al. 20031 ). as shown in Fig.EJ al- 
though in this case the implied distance from the 
JHK color-color and color-magnitude diagrams is 
106 ± 7 pc for E B -v = 0.02. 

Polaris is also argued to be a member of th e 
Pleiades moving group ( Engle fc Guinan 20041 ). 
for which the main-sequence turnoff mass is 5.9 
Mq. We have been unable to confirm a similar- 
ity in U, V, and W motions between Polaris and 
Pleiades members from an independent analysis, 
so it is unclear if such a connection exists, despite 
the close agreement in implied evolutionary ages 
for Polaris and Pleiades members. 

6. Polaris Kin? 

Stotherd (|2009h argues that Polaris and a few 
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Table 1 
Mass estimates for Polaris 



M/Mq 


Pulsation Mode 


Technique 


3.9 ±2.9 


FM 


logo and PR relation (Turner & Burke 2002) 


6.6 ±4.9 


OT 


logo and PR relation (Turner & Burke 2002) 


5.1 ±0.5 


FM 


Cluster PM relation (Turner 1996) 


6.1 ± 0.6 


OT 


Cluster PM relation (Turner 1996) 


4.5 ± 1.8 


OT? 


Hivvarcos tt and orbit (Evans et al. 2008) 
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Fig. 4. — BV photometry from Hipparcos for stars 
within 3° of Polaris. The ZAMS is plotted for d = 99 
pc and Eb-v = 0.02, along with an isochrone for t — 
8 x 10 7 yrs. 



other unusual Cepheids (e.g., V473 Lyr) may con- 
stitute a special class of "Blazhko Cepheids," on 
the basis of their unusual period and amplitude 
changes, akin to the Blazhko Effect in RR Lyrae 
variables. In this mechanism the growth and de- 
cay of a strong magnetic field modulates surface 
convection, thereby modulating pulsation ampli- 
tude and period. For stars like Polaris the pulsa- 
tion period should increase as amplitude declines 
(stronger convection), and decrease as amplitude 
rises (weaker convection). In Polaris the period 
changes are in the exact opposite sense, so the 
glitch and subsequent decline in light amplitude 
for the Cepheid over the past thirty years remain 
unexplained by such a mechanism. In fact, the 
dominant period changes observed in the large ma- 
jority of Cepheids are readily explained by the evo- 
lutionary changes in radius that occur as massive 




Fig. 5.— JHK photometry (jCutri et al] 12003! ) for 
stars in Fig. [4] The ZAMS in this case is plotted for 
d = 106 pc for E B -v = 0.02. 



stars evolve through the Cepheid instability strip 
([Turner et al.ll2006df ). 

It is possible to map the rate of period change 
diagram {P— logP) from stellar evolutionary mod- 
els for individual crossings of the Cepheid insta- 
bility strip, as in Fig. [6l The rapid first cross- 
ing occurs when the stars are undergoing rapid 
core contraction supported only by a surrounding, 
thin, hydrogen-burning shell. The crossing occurs 
at a rate almost two orders of magnitude faster 
than the slower second and third crossings, when 
the stars are undergoing core helium-burning, so a 
random sample of Cepheids should contain only a 
few objects in the first crossing, perhaps only a few 
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percent at most. Stellar evolutionary models pre- 
dict specific rates of period increase for first cross- 
ings, identified by the gray area in Fig. [6] Three 
Cepheids fall within that region: DX Gem (P = 
3 d . 14), BY Cas (P = 3 d . 22), and the double-mode 
pulsator HDE 344787 (P1/P0 = 3 d .80/5 d .40). 
DX Gem and BY Cas are suspected FM pulsators, 
the former as a possible outlying member of an 
anonymous open cluster. BY Cas lies near the 
cluster NGC 663, but seems unlikely to be a mem- 
ber. 
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Fig. 6. — Observed rates of period change for 
Cepheids (|Turner et al.ll2006dl ). displaying the likely 
association with instability crossing mode: third, fifth 
(top), second, fourth (bottom). The shaded area cor- 
responds to evolutionary model predictions for a first 
crossing. Data points corresponding to Polaris and 
HDE 344787 are identified, that for the most recent 
rate of period increase for Polaris as an open circle 
(below the filled circle datum) . Open diamonds repre- 
sent where the corresponding period changes of Polaris 
would fall if it is an OT pulsator. 

The rates of period change for DX Gem, BY 
Cas, and HDE 344787 are consistent with expecta- 
tions for a first crossing of the instability strip. Is 



that also true for Polaris? All four Cepheids have 
common characteristics: rapid period increases, 
small light amplitudes (AB = 0.49 for DX Gem, 
AB = 0.53 for BY Cas, AB = 0.06 for Polaris, 
and AB ~ 0.02 for HDE 344787), and sinusoidal 
light curves. The putative first crossers also lie 
near the middle of the instability strip, according 
to their derived reddenings (Fig. [3 That is con- 
sistent with model predictions from I Alibert et all 
( 19991 ) for stars lying towards the cool edge of 
the instability strip for first crossers. Alterna- 
tively, if the location of Polaris in Fig. [7] is ad- 
justed upwards by, say, m .6 to account for p ossi- 
ble overtone pulsation (see lTurner et al.ll2005T ). an 
anomaly arises: the star would still lie redward of 
the hot edge of the instability strip, but in a region 
where the largest amplitude pulsators are located 
(IPel fc Ludll978l : iTurner et al.ll2006dh . There is 
very little leeway possible in such an interpreta- 
tion, given the well-established spectral type and 
interstellar reddening of Polaris. 
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Fig. 7. — The location of the putative first crossing 
Cepheids in the instability strip (shaded region) as FM 
pulsat ors. The str ip boundaries are observationally 
based (|Turnerll2006h . 



The interpretation of DX Gem, BY Cas, and 
Polaris as putative first crossers is not without 
problems. DX Gem has an inferred rate of period 
increase of P = 3 S .87±0 S .19 yr _1 , consistent with 
a first crossin g, but with strange O-C wanderings 
(|Hackdll989l ) and possible amplitude changes that 
may display a slow overall decrease with time. 
Its characteristics may match those proposed by 
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Stothers for "Blazhko Cepheids" (|Stotherdl200"9l) . 
BY Cas has an inferred rate of period increase 
of P = 4 s . 60 ± s . 13 yr _1 , also consistent with a 
first crossing. But its O-C variations are in need 
of confirmation, since both recent and archival ob- 
servations display peculiarities suggesting that the 
rate of period change undergoes modulations with 
time. Possible temporal changes in its light am- 
plitude have never been investigated. In the case 
of Polaris the observed rate of period increase is 
only marginally consistent with expectations for 
a first crossing, unless the rates predicted from 
stellar evolutionary models are systematically too 
large. 

HDE 344787 is perhaps the best case for a first- 



crossing Cepheid (jTurner et al. I l2009bh . It is a 



double-mode pulsator with an observed rate of pe- 
riod increase of P — 12 s . 96 ± 2 s . 41 yr^ 1 , but also 
with a declining light amplitude (Fig. [5]) that may 
lead to a compl ete cessation of pulsa tion within a 
few years' time ( Turner et al. l2009ah . 
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Fig. 8. — The declining light amplitude of HDE 
344787 for FM puls ation (open circles) and OT pulsa- 
tion (filled circles) (|Turner et al.ll2009bl ). 

The case for OT pulsation in Polaris may also 
be addressed observationally by another method. 
There is a well-known phase lag, Atjii, between 
maximum brightness in Cepheid light curves and 
the corresponding phase of velocity minimum. 
Qgloza fe Moskalikl (|2000l ) have found that the ob- 
served period dependence of the values of A<fii 
for OT pulsators appears to differ from that of 
FM pulsators, thereby enabling one to discrimi- 
nate between FM and OT pulsation in Cepheids. 
The gener al trends appear to be confirmed in 
models by ISzabo et al. ( 2007 ). According to this 
criterion, the value o f Acfci = —0.45 ± 0.06 ob- 



changes in Polaris complicate the determination 
of A0i for the C epheid, and possibly others as 
well. The study bv lTurner et al.l (|2005l ) accounted 
for such changes directly by carefully matching 
photometric and radial velocity O-C estimates 
closely adjacent temporally, and yielded a value of 
A<£i = -0.10±0.01 (-0.383± 0.027 days), consis- 
tent w ith the FM sequen ce of lOgloza fc Moskalikl 
(|2000h . as well as that of lSachkovi (j2000f ). By tbis 
criterion, Polaris would be considered as a funda- 
mental mode pulsator. 

As is ever the case with Polaris, the situation 
is somewhat more complicated. Radial velocities 
also exist for the newly -discovered small am plitude 
Cepheid HDE 344787 (jTurner et al.ll2009rj) . which 
appears to display different phase lags for the FM 
and OT modes, although the value obtained for 
FM pulsation is uncertain because the FM signal 
is much weaker than the OT mode signal. For 
the FM pulsation in HDE 344787 (P = 5 d .40), 
Ac/)i = —0.32, whereas for the OT pulsation 
(P = 3 d .80), A<j>i = -0.09. The value of Afc for 
the fundamental mode is roughly consistent with 
the trend found bvlQgloza fc Moskalikl |7oOO), but 
not with that of Sachkovl (|2000h . whereas the value 
of A(f>i for the overton e mode disagrees markedly 
wit h the trends of bothlQgloza fc Moskalikl (|200of) 
and ISzabo et al. ( 2007 ). If the corresponding sit- 
uation for Polaris is decided in empirical fashion, 
i.e., by comparison with HDE 344787, then the 
similarity of the value of A(f>\ for OT pulsation in 
HDE 344787 with the value of Ac/>i = -0.10 for 
Polaris, which has a similar pulsation period, sug- 
gests that Polaris is likely to be an OT pulsator. 

Or is it necessarily true that double-mode 
Cepheids actually differ from single-mode pul- 
sators in their A<pi properties? It may be that 
there is no distinction in Atfii between FM and 
OT pulsation, but further careful study of the 
observational sample is needed to confirm that 
possibility. The contradictions arising from the 
arguments presented here clearly argue for the 
need of improved observational data on A^i for 
Cepheids of known pulsation mode, in particular 
data that implicitly account for the rapid period 
changes in many of the variables, a factor that 
directly affects the determination of reliable phase 
shifts for them. 



tained for Polaris by iMoskalik fc Qgloza I (|2000l ) 
implies overtone pulsation. But the rapid period 



7 



7. The Polaris Multiple System 

The recent reso l ution of the Polaris A sys- 
tem bv lEvans et al.l ( 2008 ) using the Hubble Space 
Telescope (HST) has made it possible to study 
the three previously-recognized components of the 
system directly. The two measured HST positions 
for Polaris Ab relative to the Cepheid are consis- 
ten t with the astr o metri c orbital solution proposed 
by IWielen et al. ( 2000l ). in which a retrograde 
orbit with an orbital inclination of i = 130°. 2 



was o btained. An orbital solution bv lEvans et al 



(2008) combined HST observations with astromet- 
ric data to yield the parameters in TableUJ Match- 
ing the results to stellar evolutionary models pre- 
sented some inco nsistencies w i th the inferred mass 
of the Cepheid (jEvans et al.l 120081 ) . possibly be- 



cause it was assumed that Polaris is in the third 
crossing of the instability strip, which, as noted 
earlier, conflicts with the observed rapid rate of 
period change. 



R 
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Fig. 9. — Seasonal mean radial velocities for Polaris 
from 1888 to 2005. Filled circles represent the actual 
data, while open circles and crosses denote the same 
data folded forwards and backwards in time by one 
or two orbital cycles, respectively. The upper diagram 
illustrates the confusion introduced by literature veloc- 
ities that do not appear to match the velocity system 
of Lick Observatory exactly. The lower diagram illus- 
trates the data following the adoption of corrections 
to adjust the velocities to the Lick system. 

It is unclear how the results depend upon 
the orbital solution for the sp ectroscopic subsys- 
tem. The Evans et al.l ( 20081 ) analysis of Polaris 



scopic binary from Kampei ( 19961 ). but a new 
solution for the system has been o btained using 



more r ecent radial velocity measures ([Turner et al 
2006a ). The new orbital solution depends upon 



small zero-point adjustments to existing r adial ve- 
locitie s, m ost notably the observatio ns by I Kampei 
(|l996l ) and lKamper fc Ferniel (|l998h . as illustrated 
in the upper portion of Fig. [9l Correction for 
such zero-point offsets reduces the scatter in the 
velocity observations from one cycle to another 
(Fig. [9l lower), but leaves residual scatter that 
is frequently larger than expected from the un- 
certainties in the measurements. The resulting 
phased best data set is shown in Fig. [TU] and is 
summarized in Table [7J Three separate solutions 
were made by the Lehmann-Filhes technique: one 
tied to recent observations, a second to the full 
data set (which gives a better match to the veloc- 
ity amplitude K), and a third restricted to data 
that give the best overall visual match to the 
main velocity trend. An independent solution was 
made with the same data usi ng a Bayesian Marko v 
Chain Monte Carlo analysis ([Turner et al. 2006c). 
The results, shown in the last column of Table El 
confirm the alternate solution based on the same 



data set ( Turner et al.l l2006c). 




A adopted an earlier solution for the spectro- 



0.5 1.0 
Phase 

Fig. 10. — The new orbital solution for Polaris A, 
shown as the solid curve in the diagram. Filled circles 
denote post-1970 data, open circles pre-1970 data. 

The new solution implies a less elliptical and 
larger orbit than was found in previous studies, 
but recent observations display kinks in the radial 
velocity curve for the orbit near velocity maximum 
that argue for the need of additional velocity cov- 
erage. It is possible to use the new orbital solution 
with the recent visual resolution of the compan- 
ion ([Evans et al.l 12008) to constrain the distance 
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Table 2 

Orbital Parameters for Polaris A (|Evans et al.I|2008T ) 



Parameter Value 

Orbital inclination, i 128° ± 21° 

Ascending node, fi 19° ± 15° 

Semi-major axis, a 0".133 ± 0".015 
Total system mass, 

Primary mass, M Aa ^-i'a M @ 

Secondary mass, M&b 1.26to'_7 Mq 



to the system. For this analysis the mass of Po- 
laris is assumed to lie between 3 and 6 Mq, while 
the mass of the F6 V companion is assu med to 



be less than 1.45 M (jEvans et alJl2008h . Such 



constraints produce the results shown in Fig. [TT1 
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Fig. 11. — Permitted parameters for the Polaris sys- 
tem (the unshaded region) according to the orbital 
solution and HST detection of the companion. The 
sloped curves denote, from top to bottom in each plot, 
orbital inclinations of 50°, 60°, 70°, and 90°, respec- 
tively. Shaded regions denote implausible orbital solu- 
tions dictated by the observed faintness of the compan- 
ion (upper cutoff) and the limiting orbital inclination 
of i = 90° (lower cutoff). 

The results of the analysis can be summarized 
as follows. The allowable range of parameters 
yields (d) = 121 ±10 pc and (M v ) = -3.49 ±0.18 
for Polaris, consistent with the Hipparcos paral- 



lax ( van Leeuwen et al. 20071 ) and overtone pul- 
sation. However, the most likely orbital inclina- 
tion is i > 60°, which conflic ts with the astro - 
metric estimate of i = 50° ( Wielen et al. 2000) 
unless Mp i ar i s < 3.5M . The permitted so- 
lutions are inconsistent with the inferred dis- 



tance of the visual companion (Polaris B ) (jTurner 
1977L l2006t lUsenko fc Klochkoval 120081) . the ob- 



served rapid rate o f evolutionary period change 
(jTurner et al.ll2005l ). the lack of any signal tied to 
putative FM pulsation around 5 d . 7 in the light 



and radial velocity observ a tions (iTurner 
Sprecklev fc Stevens! |2008t iBruntt et all 



2006; 



2008; 



Lee et al. 20081) . and possible me mbership in th e 



surrounding anonymous cluster (jTurner 120061) . 
The preferred solutions cluster near i = 90°, al- 
though an edge-on orbit is inconsiste nt with the 
HST observations (jEvans et al.ll2008l) . In other 



words, no solution for the system parameters is 
ideal. 

Another approach can be made using the O-C 
data of Fig. __J since deviations from the O-C fit 
should also yield parameters for the binary system. 
The light travel time delays for a star orbiting in 
a binary system depend upon a\ sin i in similar 
fashion to the radial velocities, although they are 
a quarter cycle out of phase. The largest weight 
O-C observations for Polaris do indeed follow such 
a trend, as shown in Fig. []__] What is remarkable 
is the magnitude of the deviations, which seem 
to reach almost 0.15 light day, or 26 A.U., a value 
that seems inordinately large for the motion of the 
most massive star in the system about its (presum- 
ably) less massive companion. 

The last result can be seen in greater detail in 
Fig. [T__| which plots the O-C residuals as a func- 
tion of orbital phase relative to expectations from 
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Table 3 

Orbital Solutions for Polaris A 



Parameter 


Roemer (1965) 


Kamper (1996) 


Turner et al. (2006c) 


Turner et al. (2006e) 








BMCMC Analysis 


Period, P (years) 


30.46 ± 0.10 


29.59 ± 0.02 


29.71 ± 0.09 


29.80 ± 0.05 


(days) 






10852 ± 33 


10885 ± 17 


Orbital eccentricity, e 


0.639 ± 0.012 


0.608 ± 0.005 


0.543 ± 0.010 


0.52 ± 0.01 


Epoch of periastron 


1928.48 ± 0.08 


1928.48 ± 0.08 


1928.57 ± 0.06 




JD 






2425453 ± 22 


2425430 ± 44 


Longitude of node, u? (°) 


307.24 ± 1.82 


303.01 ± 0.75 


309.6 ± 0.7 


301 ± 2 


Semi-amplitude, K (km s -1 ) 


4.09 ± 0.10 


3.72 ± 0.03 


4.41 ± 0.07 


4.23 ± 0.07 


Systemic velocity, 7 (km s~ 1 ) 


-16.41 


-16.42 ± 0.03 


-15.90 ± 0.06 


-15.87 ± 0.05 


Projected semi-major axis, ai sin i (A.U.) 


3.22 ± 0.10 


2.90 ±0.03 


3.69 ± 0.09 


3.62 ± 0.06 



A(O-C) 
(days) 



+ 5 







a UMi 





241 242 243 244 245 246 
HJD/ 10000 

Fig. 12. — The deviations of the high quality O-C 
data from the best-fitting parabolic trends matched 
to the observations. Tick marks denote epochs when 
the Cepheid primary should be crossing into the plane 
of the sky. 



the m ost recent orbital solution ((Turner et al 
2006c) . It should be noted that some of the most 



deviant trends in Fig. [13] are tied to recent obser- 
vations of the Cepheid from space-based platforms 



(|Sprecklev fc Stevend l2008: iBruntt et al . 2008), so 
it is difficult to attribute them to "observational 
error." A similar problem has been detected in 
O-C observations o f the binary Cepheid RT Aur 
([Turner et al.l l2007h . although even in that case 
there is no satisfactory resolution of the problem 
by tying it to problems with the observations. Pos- 
sibly the effect is simply another manifestation of 
the random fluctuations in pulsatio n period that 



are de tected in almost all Cepheids ([Turner et al 



2009d) , although the trends in the O-C deviations 
seem to be too well organized for a "random" 
process. 

8. Additional Considerations 

The numerous contradictions noted here re- 
garding the pulsation mode of Polaris suggest the 
possibility of a more fundamental problem. Some 
of the anomalies in the radial velocity measures 



0.5 
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(days) 
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Fig. 13. — Phased deviations of the O-C data from 
the best-fitting parabolic trends, shown relative to 
the orbital solution of Table [7] (solid line), with the 
size of the symbol increasing in proportion to its as- 
signed weight. A dashed line denotes expectations for 
ai sini = 5x the actual value. 



may be a consequence of different sets of spectral 
lines (ionization and excitation states) and wave- 
length coverage used to infer the photospheric mo- 
tion, or to the inevitable problems of establishing 
reliable systemic velocities for a pulsating star, 
where shock effects are d etected at some phases 
( Hatzes fc Cochran l2000h . But that does not re- 
solve the problems arising from the O-C offsets 
for the photometric observations, which are of rel- 
atively high precision for recent data sets. The 
presence of a third star in the spectroscopic sub- 
system, as surmised by the re sults of most previous 



radial velocity studies (e .g., Dinsha w et al. 1989; 



Hatzes fc Cochran! l2000l: iLee et al.ll2008l ). might 



resolve the various discrepancies in our current un- 
derstanding of Polaris as a Cepheid, but would 
require continued spectroscopic and astrometric 
monitoring of the star for some years to come in 
order to provide confirmation. It is interesting to 
note that the "glitch" in the O-C observations 
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for Polaris occurred only a few years after peri- 
astron passage for the spectroscopic subsystem in 
1958, when the stars were near closest approach 
to one another. Presumably that would be the 
most likely interval of time when an interaction 
with another close companion would occur. But, 
of course, that is merely speculation. 

The reasonable possibility that some of the 
stars near Polaris might be ph ysically assoc iated 
with the Cepheid was noted by Ferniei (1966), but 
has not been pursued to a great extent, possi- 
bly because of the difficulties of observing stars 
in close proximity to the north celestial pole us- 
ing ground-based facilities. The lack of detectable 
X-ray emission from the faint red stars C and D 
in the Polar i s mul tiple system has been used by 
Evans et all (|2006h as evidence that they are not 
young, newly- arrived, main sequence dwarfs, as 
required by a possible association with Polaris, 
so are most likely field stars. There are several 
other brighter, potential members of the Polaris 
cluster lying outside the immediate field of the 
Cepheid, however, for which basic photometric 
and radial velocity observations do not exist. A 
program to obtain such data has been underway 
for a few years now, but the study is far from com- 
plete. Other low-amplitude sinusoidal Cepheids, 
the "Polaris kin," are even less well-studied, al- 
though they could prove to be equally valuable 
for clarifying the nature of Polaris itself. Clearly 
the Polaris story is not yet over. 
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